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<D . Abstract 

We re-analyse the resonant spin-flavour (RSF) solutions to the solar neutrino 
^ ' problem in the framework of analytic solutions to the solar magneto-hydrodynamics 

h : 

(MHD) equations. By substantially eliminating the arbitrariness associated to 
the magnetic field profile due to both mathematical consistency and physical re- 
quirements we propose the simplest scheme (MHD-RSF, for short) for solar neu- 
trino conversion using realistic static MHD solutions. Using such effective two- 
parameter scheme we perform the first global fit of the recent solar neutrino data, 
including event rates as well as zenith angle distributions and recoil electron spectra 
induced by solar neutrino interactions in Superkamiokande. We compare quantita- 
tively our simplest MHD-RSF fit with vacuum oscillation (VAC) and MSW-type 
(SMA, LMA and LOW) solutions to the solar neutrino problem using a common 



1 E-mail: omr@fis.cinvestav.mx 

2 E-mail: penya@flamenco.ific.uv.es 

3 E-mail: rashba@izmiran.rssi.ru 

4 E-mail: semikoz@flamenco.ific.uv.es 

5 E-mail: valle@flamenco.ific.uv.es 



well-calibrated theoretical calculation and fit procedure. We find our MHD-RSF 
fit to be somewhat better than those obtained for the favored neutrino oscillation 
solutions, though not in a statistically significant way with Am 2 « 10~ 8 eV 2 and 
sin 2 26 = 0. We briefly discuss the prospects to disentangle our MHD-RSF sce- 
nario from oscillation-type solutions to the solar neutrino problem at future solar 
neutrino experiments, giving some predictions for the SNO experiment. 
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1 Introduction 



The persistent disagreement between solar neutrino data and theoretical expectations 
has been a long-standing problem in physics. Since the very first measurements the 
Solar neutrino problem has remained as a puzzle, re-confirmed by new data on rates by 
GALLEX-SAGE |]2], [| as well as most recently published 825-day data collected by the 
Super-Kamiokande collaboration Q which goes beyond the simple rate measurements to 
include also rate-independent data such as the recoil electron spectra induced by solar 
neutrino interactions, as well as the zenith angle distributions J5J. It has often been 
argued that these data can not be accounted for by astrophysics Together with the 
atmospheric neutrino data |7| these constitute the only present-day evidence in favour 
of physics beyond the Standard Model, providing a strong hint for neutrino conversion. 

The most popular solutions of the solar neutrino anomalies are based on the idea 
of neutrino oscillations, either in vacuum or in the Sun due to the enhancement arising 
from matter effects ||. 

Although these are the simplest neutrino conversion mechanisms there is consider- 
able interest in alternative interpretations. For example it has long been noted |§ that 
Majorana neutrinos may have non-zero transition magnetic moments which can gener- 
ate spin-flavour conversions in the presence of a magnetic field. These are especially 



interesting for two reasons: (i) on general grounds [ 10j neutrinos are expected to be 
Majorana particles and (ii) conversions induced by transition magnetic moments can be 
resonant in the Sun fllTl . There is also room for more exotic mechanisms such as flavour 



changing neutrino interactions |12j which do not even require neutrino mass fl"3|j . 

Here we will re-analyse the status of resonant spin-flavour solutions to the solar 
neutrino problem in the light of the most recent global set of solar neutrino data, 
including event rates as well as zenith angle distributions and recoil electron spectra 
induced by solar neutrino interactions in Superkamiokande which has attracted inter- 
est recently |T3|, [13, HJ. In contrast to previous attempts we will adopt the general 



framework of self-consistent magneto-hydrodynamic (MHD) models of the Sun fl7 |. 
A previous attempt in this direction is given in ref. |18|]. For definiteness we will con- 
centrate in the recent proposal of Ref. |TP[ where relatively simple analytic solutions 
have been given. We perform global fits of solar neutrino data for realistic solutions 
to the magneto-hydrodynamics equations inside the Sun. This requires adjusting both 
the neutrino parameters as well as optimizing the magnetic field profile. The arbitrari- 
ness associated to the latter is substantially reduced due to mathematics (they must 
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be solutions of MHD equations) as well as reasonable physical requirements. This way 
and by neglecting neutrino mixing we obtain the simplest MHD-RSF solution to the 
solar neutrino problem, characterized by two effective parameters, Am 2 and [i v B± max , 
B±max being the maximum magnitude of the magnetic field inside the convective region. 
Throughout this paper we have assumed that the neutrino transition magnetic moment 
\i v is given in units of fin = fi u /10~ u /j,b, where \ib is the Bohr magneton and we set 
/in = 1 everywhere. Our MHD-RSF solution can be meaningfully compared with the 
neutrino oscillation solutions to the solar neutrino problem. We find that our simplest 
two-parameter MHD-RSF fits to the solar neutrino data are slightly better than those for 
the oscillation solutions, but not in a statistically significant way. The required best fit 
points correspond to maximum magnetic field magnitudes in the convective zone smaller 
than 100 KG. We briefly discuss the prospects to distinguish our simplest MHD-RSF 
scenario from the neutrino oscillation solutions to the solar neutrino problem at future 
solar neutrino experiments, giving some predictions for the SNO experiment. 



2 Static Magnetic Field Profiles in the Sun 



In solar magneto-hydrodynamics |3(J (MHD, for short) one can explain the origin of 
solar magnetic fields from the dynamo mechanism at the bottom of the convective zone 
or, to be more specific, in the overshoot layer, where magnetic fields may be as strong 



as 300 kG ||21|| . Such a picture is quite attractive and several MHD dynamo solutions 



has been known since long time ago (see for example [T7|]) However the corresponding 
magnetic field profiles are rather complicated and difficult to extract. For this reason 
there have been many attempts to mimic MHD properties through the use of ad hoc 
magnetic field profiles involving, for example, twisting fields [[22]. 

Here we will follow an alternative approach using fully self-consistent solutions to the 
MHD equations inside the Sun. To achieve this we focus on the case of stationary solu- 



tions which are known analytically in terms of relatively simple functions fL9 |. This way 
we obtain a simple and well-motivated magnetic field profile, without the full complexity 
that a dynamo model implies. In this section we will explain this model and discuss the 
limits on the shape parameters describing the field profile. We will also discuss how to 
relate this model with the dynamo picture of the solar interior. 
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2.1 Single-Mode Field Configurations 

In this subsection we will describe the model we are using for the magnetic field profile. 
We consider only solutions to the equation for a static MHD plasma configuration in a 
gravitational field given by the equilibrium of the pressure force, the Lorentz force and 
the gravitational force 

Vp- -j x B + pV<£> = 0, (1) 

where p is the pressure, j = (c/47r)roi B is the electric current, B is the static magnetic 
field under consideration, p is the matter density and $ is the gravitational potential. 

This static MHD equations correspond to a quiet Sun and they admit axially sym- 
metric solutions in the spherically symmetric gravitational field which can be simply 
expressed in terms of spherical Bessel functions and were first discussed in Ref. |H| . 
For this model the magnetic field will be given by a family of solutions that depends on 
Zk, the roots of the spherical Bessel function f 5 / 2 = \fzJ$n{z)-> to ensure the boundary 
condition that B vanishes on the solar surface. Within the solar interior the magnetic 
field for any k will be then given by 
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Here 9 is the polar angle and the distance r has been normalized to R & = 1. Taking into 
account the inclination of the solar equator to the ecliptics, where neutrinos propagate to 
the Earth, it follows that 9 lies in the narrow range 83° — 97°, depending on the season. 
In our calculations we have averaged over 9 in the above range. 

The modulus of the perpendicular component which is relevant to the neutrino spin- 
flavour takes the form 

sin^ 



Si 



Bl + Bn = B c 



-f(r) , 



(4) 



where f(r) is some known smooth function. Notice also that the behaviour of B at the 
solar center (r = 0) 
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Figure 1: The perpendicular component of B for various A;- values 1 (solid), 3 (dashed) 
and 10 (dotted). 

B g (0,e) = -B core sm6 , 

Bf(0,e) = 5 core sin^-^->0, (5) 

is completely regular, determined only by the parameter B core . In Fig. 1 we display the 
perpendicular component of B for various fc-values 1, 3 and 10, which correspond to the 
roots z k = 5.7, Zk = 12.3 and z k = 34.5, respectively. 

2.2 Astrophysical Constraints on Magnetic Fields 

We now discuss the astrophysical restrictions on the free parameters B core and k char- 
acterizing the model. We can see that the magnitude of a magnetic field at the center 



of the Sun is constrained by the Fermi- Chandrasekhar limit [23| which implies [19 



2 5 7 Q-6 {B k fzl k R% 
15 to - 1 GM? 



< 1, (6) 



'0 

here 70 is polytropic index characterizing the equation of state (pressure P ~ p 10 )-, M & 
is the solar mass and G is Newton's constant. This equation gives us an upper bound 
on B core < 2 MGauss for k = 1. For higher values of k this constraint is even weaker. 

Regarding with the values of k. These can be constrained by taking into account 
that in order to justify the use of a stationary solution, it is necessary that the diffusion 
time due to ohmic dissipation 

4ncr cond L 2 (k) 

tdiss{k) = 5 (7) 



must be bigger than the age of the Sun t & ~ 1.4 x 10 17 s |33| . Here L(k) denotes a 
characteristic spatial scale of the magnetic field which corresponds to the typical distance 
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between subsequent nodes of the corresponding Bessel function. As we can see from 
Fig. 1 L(k = 1) ~ R®, while L(k) ~ R Q /k. In eq. ([?]) <r cond = uj 2 / \A.iiv ep ) denotes 
the conductivity of the fully ionized hydrogen plasma. After substituting the plasma 
frequency 



w„ = 5.65 x 10 4 ,/ Ue - s- 1 



and the e-p collision frequency 



lcm -3 



Ve V = 50 ( ] I — ] s 



,lcm~ 3 / \1K J 

we obtain from eq. (0) an estimate of the magnetic field dissipation time tdi SS {k) 

t d iss(k) = 1 x 1 '' > t Q (8) 

c z k z 

Note that the dissipation time is shorter for higher k values, long-lived field configura- 
tions being possible only if k is small. For example, the dissipation time for the third 
mode is about an order magnitude less than that for the first mode. Moreover it depends 
on the value of the temperature that we take. Some typical values for the temperature 
are T mm ~ 2.8 x 10 5 if for the bottom of the convective zone and T max ~ 1.6 x 10 7 K for 



the solar core p5| . Thus, taking the optimistic estimate, T max , we obtain k < ku = 13, 
while, if we consider the average value T = 4 x 10 6 K we will have k < k]M — 5. In what 
follows we will consider values of k < 10. 



2.3 Energy localization criterium 

It is commonly accepted that magnetic fields measured at the surface of the Sun are 
weaker than within the convective zone interior where this field is supposed to be gen- 
erated. It is known by observational data that the mean field value over the solar disk 
is of the order of 1 Gauss while in the solar spots magnetic field strength reaches 1 KG. 

On the other hand the general knowledge of the solar magnetic field models is that 
the magnetic field increases at the overshoot layer, while being small at the solar interior, 
a picture rather opposite to the one we have seen in Fig. 1. 

Although there is no direct information on the magnetic field magnitudes at the 
solar core, there are theoretical reasons which imply a central magnetic field less than 



30 Gauss f26j; otherwise the present magnetic field in the convective zone would be too 
big, leading to a visible enhancement in sunspot activity. 

This conflict can be avoided by taking advantage of the linear nature of the basic 
equilibrium MHD equation in eq. (|1|). This implies that any linear combination of 
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Figure 2: Magnetic field configurations obtained by combining individual modes for 
different ku values, 5, 6 and 10. Summing up to higher modes achieves better localization 
of the field in the convective region (solid). 

solutions B k {k = 1, 2, ... , k^, for some fixed number kja) 

B = ci-Bi + c 2 B 2 + ... + c M B M (9) 

is also a solution. As mentioned in section 2.2 we will adopt ku < 10 in order to ensure 
that ohmic dissipation is acceptable and therefore justify the static approximation. 

In order to ensure that the magnetic field energy is localized mainly within the 
convective region we will now supplement the constraints of section 2.2 by imposing that 
the magnetic field should vanish in the center of the Sun [] 

B(f=0)=0 (10) 

The latter implies 

C! + c 2 + ... + c kM =0. (11) 

Therefore we will have, in principle, ku — 1 free parameters. 

We will require, in addition, that the magnetic field energy must be minimal in the 
region below the bottom of the convective zone, characterized by a certain value of r$ 

JO 07T 

This implies 

-^E B = £)ff Vii = (13) 



6 Note that adopting a finite but small value for \B(f = 0)| < 30 Gauss does not change significantly 
the profile of the total magnetic field eq. (^) which is the relevant quantity for describing neutrino 
propagation. 
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where 



aij = J d 3 r(Bi - B M ) ■ B* . (14) 



Without loss of generality we can assume that one of the coefficients is non-zero, 
which prevents us from having only the trivial solution q = 0. Taking Ck M -i 7^ we will 
have the linear non-homogeneous equation 

Zf^Cjdij = -c kM . 1 a ikM _ 1 (15) 

which determines all the q coefficients in terms of, say, Ck M -i- As expected on physical 
grounds, this last remaining parameter Ck M -i corresponds to the maximum magnetic 
field magnitude in the convective region, i.e. B± max is proportional to Ck M -i- In general, 
all the coeficients Cj will be different from zero for i < M and alternate in sign. As an 
example for k M = 6 the coefficients are -0.338968, -0.261825, 1.29186, -1.77360, 1.86916, 
-0.786628. 

The procedure sketched above provides a consistent method for combining individual 
mode solutions By. of the static MHD equation, while fixing all of the coefficients of 
the linear combination, leaving as free parameters only the value of B± max inside the 
convective region and the value of k^i — 10. In Fig. (0) we show the resulting combined 
profiles for ku = 5, 6, 10. The parameter r could also be taken as a free parameter but, 
on physical grounds, it should lie in a narrow range close to overshoot layer. We show 
explicitly that varying r has little effect on our results. 



3 Fitting the Solar Neutrino Data 

The Majorana neutrino evolution Hamiltonian in a magnetic field is well-known to be 
four-dimensional |5J. For definiteness and simplicity we will neglect neutrino mixing in 
what follows and consider first the case of active-active neutrino conversions. This will 



allow us to compare our x -analysis with the previous ones [14, 15, 16]. The v e — > vi 






-( 


U / 





conversions are described by the master Schrodinger evolution equation 

V e — 5 u v B + \ ( v e \ 

(16) 

\l v B- -Vi + 5 J y P £ J 

where \x v denotes the neutrino transition magnetic moment in units of 10~ n £ 
denoting either fi or r. Here B± = B x ± iB y and 5 = Am 2 /AE is the neutrino mass 
parameter; V e (t) = G F V2(p(t)/m p )(Y e - Y n /2) and V e {t) = G F V2{p(t)/m p ){-Y n /2) 
are the neutrino vector potentials for v e and vi in the Sun given by the abundances of 
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Figure 3: Typical MHD-RSF neutrino survival probablity P ee versus E/Am 2 . 

the electron (Y e = m p N e (t) / p(t)) and neutron (Y n = m p N n (t) / p{t)) components. In 
our numerical study of solar neutrino data we adopt the Standard Solar Model density 
profile of ref . . 

We solve Eq. (jTGj) numerically by finding a solution of the Cauchy problem in the 
form of a set of wave functions u a (t) =\ v a (t) | e** a ^ from which the neutrino survival 
probabilities P aa if) = v*v a are calculated. They obey the unitarity condition ^ a Paa = 1 
where the subscript a denotes a = e for v e and a = £ for \>i respectively. 

As an illustration we display in Fig. (|3|) the electron neutrino survival probablity 
P ee calculated in the MHD-RSF scheme from eq. (p~6f ) plotted versus E/Am 2 . This is 
obtained with the magnetic field configurations given in Fig. (0). 

We will now re-analyse the status of resonant spin-flavour solutions to the solar 
neutrino problem in the light of the most recent global set of solar neutrino data, in- 
cluding event rates as well as zenith angle distributions and recoil electron spectra in- 
duced by solar neutrino interactions in Superkamiokande which has attracted interest 
recently [[U], |l^, |l6l . It has been found that the quality of the fit to the solar neutrino 



data depends on the magnetic field profile. The best solutions have been obtained with 
a magnetic field around 100 KG in the convective zone and zero at the core (profile 3 
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Experiment 


Rate 


Ref. 


Units 


^>BP98 


Homestake 
GALLEX + SAGE 
S up er-Kamiokande 


2.56 ±0.23 
72.3 ±5.6 
2.45 ± 0.08 


1 

HI 

E 


SNU 
SNU 

10 6 cm"2 s" 1 


7.8 ± 1.1 
130 ±7 
5.2 ±0.9 



Table 1: Solar neutrino rates measured in the Chlorine, Gallium and Super-Kamiokande 
experiments. 



in Ref [[14J], and profile 6 of Ref [|15|]) or an almost constant magnitude of the magnetic 



field, but with twisting direction ||T6[ using the profiles given in |[22j| . 

In contrast with previous work we will consider the fits obtained when we employ self- 
consistent solutions of MHD equations which obey the physical requirements we derived 
in sec. 2.2 using the procedure for combining magnetic field modes described in section 
3. Our approach is global and allows us to compare quantitatively with other solutions 
to the solar neutrino problem within the same well calibrated theoretical calculation and 
fit procedure. 



3.1 Rates 

In order to determine the possible values of the parameters characterizing the MHD-RSF 
solution to the solar neutrino problem, we have first used the data on the total event rates 
measured at the Chlorine experiment in Homestake |I[ , at the two Gallium experiments 
GALLEX and SAGE || |3j and the 825-day Super-Kamiokande data sample, as given 
in table [|. 

In our statistical treatment of the data for the combined fit we adopt the \ 2 definition 



given in ref. [ Z7J , 

xl = E {Rf - RT P ) , rf(R? - Rf p ) (17) 

i,j=l,3 

where Rj h is the theoretical prediction of the event rate in detector i and Rf p is the 
measured rate. The error matrix contains not only the theoretical uncertainties but 
also the experimental errors, both systematic and statistical. 

The general expression of the expected event rate in the presence of oscillations in 
experiment % is given by Rf 1 



Rf = E fa fdE u X k (E u ) x [a e>i {E v ) (P ee (E u ,t)) 

fe=l,8 

+<r x ,i(E u )(l- (P ee (E v ,t)))}, 
10 




Figure 4: MHD-RSF 90% CL (light) and 99% CL (dark) regions of Am 2 versus 
B± m ax(KG) allowed by the rates given in table p], for r = 0.6 and K M = 6 

where E v is the neutrino energy, <pk is the total neutrino flux and is the neutrino 



energy spectrum (normalized to 1) from the solar nuclear reaction k [28] with the 
normalization given in Ref. |25| . Here a e ^ (cr x ^) is the v e {y x ) (with x being Ji or s 
corresponding to active-active or active-sterile MHD-RSF conversions) cross section in 



the Standard Model [^| with the target corresponding to experiment i, and (P ee (E u ,t)} 
is the time-averaged v e survival probability. 

For the Chlorine and Gallium experiments we use improved cross sections <Ji(E) from 
Ref. PJJ. For the Super-Kamiokande experiment we calculate the expected signal with 
the corrected cross section given in the Appendix Sec. |A]. 

The expected signal in the absence of oscillations, -Rf P98 , can be obtained from 
Eq. flilf ) by substituting P ee = 1. In table |l] we also give the expected rates at the 



different experiments which we obtain using the fluxes of Ref. |25| . 

In Fig. (|4]) we display the region of MHD-RSF parameters allowed by the solar 
neutrino rates. 

Our x 2 analysis of the solar neutrino rates uses the magnetic field profiles discussed 
in Section 2. As we mentioned in that section, these profiles are characterized by ku 
and ro- In table we present the best fit points for ku from 4 to 8 and for r$ ~ .6RQ 
and for Bj_ max < 300 KGauss. In the same table we also show the best fit points for 
B± max < 100 KGauss. We can see from this table that the \ 2 is pretty stable and does 
not depend significantly on the choice of ku and r allowed by astrophysics. In Fig. 
(f|) we display the region of MHD-RSF parameters allowed by the solar neutrino rates 
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0.62^0 






0.64^0 




M 




Am 2 


v 2 . 

/V TYltTl 




Am 2 


v 2 . 




Am 2 


Y 2 ■ 


4 


29. 


.84 x 10~ 8 


1.3 


32 


.94 x 10~ 8 


1.1 


41 


.99 x 10~ 8 


1.3 




29. 


.84 x 10~ 8 


1.3 


32 


.94 x 10~ 8 


1.1 


41 


.99 x 10~ 8 


1.3 


5 


71 


1.5 x 10~ 8 


1.0 


63 


1.5 x 10~ 8 


.87 


72 


1.6 x 10~ 8 


.95 




71 


1.5 x 10~ 8 


1.0 


63 


1.5 x 10~ 8 


.87 


72 


1.6 x 10~ 8 


.95 


6 


244 


1.9 x lO" 8 


.03 


247 


1.9 x 10~ 8 


.06 


251 


1.9 x lO" 8 


.13 




77 


1.1 x 10~ 8 


.33 


80 


1.0 x 10~ 8 


.40 


80 


1.0 x 10~ 8 


.35 


7 


208 


1.2 x 10~ 8 


.20 


210 


1.0 x 10~ 8 


.13 


215 


.94 x 10~ 8 


.26 




83 


.75 x 10~ 8 


.52 


84 


.71 x 10~ 8 


.54 


84 


.64 x 10~ 8 


.47 


8 


220 


.98 x 10~ 8 


.38 


222 


.94 x 10~ 8 


.24 


225 


.84 x 10~ 8 


.45 




87 


.64 x 10~ 8 


.68 


87 


.59 x 10~ 8 


.64 


87 


.55 x 10" 8 


.69 



Table 2: Best fit points for the rates-only analysis for different r and k M values in 
active-active MHD-RSF oscillations. 



for the case M = 6 and r = .QRQ. We can see that there are several allowed regions 
for different values of the magnetic field. As we already mentioned, in our analysis we 
have fixed the value of \i v to be 10 _11 /i£. Since the evolution equation depends on the 
product fi u B, for a smaller value of the neutrino magnetic moment the B^ max in Fig. 
(^D would have to be correspondingly increased. In this sense, the local minima shown 
in table for B± max < 100 KGauss allows a smaller /v 



3.2 Zenith and Spectrum Fit 

Apart from total event rates the water Cerenkov experiment also measures the zenith 
angle distribution of solar neutrino events as well as their electron recoil energy spectrum 
with their recent 825-day data sample [[|. 

The smallness of the Am 2 values indicated by the rates fit implies that no appreciable 
day-night variation of the counting rates is expected in our MHD-RSF solution. However 
the measured solar neutrino zenith angle data must be included in the analysis and we 
do that. This is necessary in order to enable us a meaningful comparison with vacuum 
and matter oscillations using the same statistical criteria |27, |3l], |32]], see definitions in 
the appendix. We obtain xlenith = 5-4 for the full range of parameters in the analysis, 
the same as for the no-oscillation case. 
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Figure 5: MHD-RSF 99% CL regions of Am 2 versus B± max forbidden by the recoil 
electron spectrum data given in table |7] of the appendix, for r = 0.6 and ku = 6 

The recoil electron energy spectrum induced by solar neutrino interactions after 504 
days of operation is given for energies above 6.5 MeV using the Low Energy (LE) analysis 
in which the recoil energy spectrum is divided into 16 bins, 15 bins of 0.5 MeV energy 
width and the last bin containing all events with energy in the range 14 MeV to 20 MeV. 
Below 6.5 MeV the background of the LE analysis increases very fast as the energy 
decreases. Super-Kamiokande has designed a new Super Low Energy (SLE) analysis in 
order to reject this background more efficiently so as to be able to lower their threshold 
down to 5.5 MeV. In their 825-day data |4| they have used the SLE method and they 
present results for two additional bins with energies between 5.5 MeV and 6.5 MeV. In 
the appendix we present these data in table [7] as well as the details of our statistical 
analysis. Our results are almost independent of the choice of the parameters ku and r 
in the physical range, thus establishing the robustness of the fit procedure. The predicted 
spectrum is essentially flat except for the upper part of the Am 2 region. As an example, 
we show in fig. |5| the excluded region at 99 % CL for the case k M = 6 and r = 0.6. 

3.3 Global Fit 

As we have seen in the partial analysis, zenith and spectrum are essentially flat in the 
region of parameters which provide a good fit for the rates-only analysis. For this 
reason, the allowed regions are slightly modified by the inclusion of the zenith angular 
dependence and the energy spectrum data. As the results are statistically independent 
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M 


B±max 


Am 2 


v 2 . 




Bj_ m ax 


Am 2 


v 2 . 




5 


72 


1.6 x 10~ 8 


25.7 


1.10 


72. 


2.0 x 10~ 8 


25.7 


1.17 




72 


1.6 x 10" 8 


25.7 


1.10 


72. 


2.0 x 10" 8 


25.7 


1.17 


6 


240 


1.8 x 10" 8 


24.7 


0.97 


241 


1.8 x 10~ 8 


24.9 


0.96 




80 


1.1 x lO" 8 


25.7 


1.11 


80 


1.1 x lO" 8 


25.9 


1.10 



Table 3: Global best fit points and local minima for B± max < 100 kGauss for different 
ro and ku values in active-active MHD-RSF conversion scenario. 



of the choice of ku and r in the physical range, our analysis effectively involves only two 
parameters. It is therefore meaningful to compare it with the popular two-neutrino fits 
characterizing vacuum or matter-enhanced oscillations. In table we show the best-fit 
points in the range of our study for different ku and r values. Moreover, we show the 
local (global) minimum for B± max less than 100 KG, which will be important to improve 
sensitivity on the transition magnetic moment of the neutrino. In fig. |5| we show the 
allowed region at 90% CL and 99% CL for the case r = 0.6 and ku = 6. We have also 
investigated the effect of varying the hep flux, obtaining for the allowed regions results 
similar to the no-oscillation solution discussed previously in ref. P?| , independently of 



the Am 2 and B± max value, with a hep normalization factor of 13.5. 

We now move to the case of active-sterile MHD-RSF conversions. For this case one 
must make substitute v s for vi in eq. ( |T6"D and take into account that V s — 0. The results 
given above for active-active MHD-RSF conversions change when conversions involve 
sterile neutrinos. The best fit points ( and local ones ) are obtained with parameters 
slightly modified with respect to those obtained for the active-active case. In the rates 
only fit, the xlates ls worse than for the active- active case, essentially due to the neutral 
current contribution in the Super-Kamiokande experiment. The zenith angle dependence 
and the recoil energy spectrum remains flat as before. The global fit for different ro and 
ku-, is shown in table §. 

4 MHD-RSF versus Oscillation Solutions 
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Figure 6: 90% CL (light) and 99% CL (dark) allowed MHD-RSF regions in Am 2 and 
B±max from the measurements of rates combined with the zenith angle distribution and 
the recoil energy spectrum in super-Kamiokande, for ro = 0.6 and ku = 6. 







O.6i? 








O.Q2R G 






M 


B± m ax 


Am 2 


y 2 . 


Ol-sp 




Am 2 


v 2 . 
Amm 




5 


68. 


1.4 x 10~ 8 


34.0 


1.25 


68 


1.5 x 10" 


-8 


35.2 


1.30 




68. 


1.4 x 10" 8 


34.0 


1.25 


68 


1.5 x 10" 


-8 


35.2 


1.30 


6 


247. 


1.9 x 10^ 8 


28.7 


1.20 


247. 


1.9 x 10" 


-8 


28.2 


1.18 




75. 


1.1 x 10~ 8 


30.5 


1.13 


75. 


1.0 x 10- 


-8 


30.6 


1.10 



Table 4: Best fit points and local minima for B± max < lOOkG for the global analysis for 
different ro and ku values in active-sterile MHD-RSF conversions. 



15 



4.1 Present 



From the results of the previous section it follows that our MHD-RSF solution to the 
solar neutrino problem provides a good description of the most recent solar neutrino data, 
including event rates as well as zenith angle distributions and recoil electron spectra 
induced by solar neutrino interactions in Superkamiokande. We have shown that our 
procedure is quite robust in the sense that the magnetic field profile has been determined 
in an essentially unique way. This effectively substitutes the neutrino mixing which 
characterizes the oscillation solutions by a single parameter B± max characterizing the 
maximum magnitude of the magnetic field inside the convective region. The value of 
kM characterizing the maximum number of individual modes superimposed in order to 
obtain a realistic profile and the parameter ro characterizing the location of the convective 
region are severely restricted. The allowed &m values are restricted by ohmic dissipation 
arguments to be lower than 10 or so, while r is close to 0.6i? Q . We have found that our 
solar neutrino fits are pretty stable as long as ku exceeds 5 and r lies in the relevant 
narrow range (see tables |^ and ^). Therefore our fits are effectively two-parameter fits 
(Am 2 and B± max ) whose quality can be meaningfully compared with that of the fits 
obtained for the favored neutrino oscillation solutions to the solar neutrino problem. In 
table |5] we compare the various solutions of the solar neutrino problem with the MHD- 
RSF solutions for the lower magnetic field presented here. Clearly the MHD-RSF fits 
seem somewhat better (though not in a statistically significant way) than those obtained 
for the MSW effect |27| as well as just-so solutions f32"f . Notice that in table [| we have 
used the same common calibrated theoretical procedures and statistical criteria. These 
results are for the case where the BP98 Standard Solar Model is adopted. We have also 
investigated the effect of varying the hep flux, obtaining a hep normalization factor of 
13.5 to be compared with 12 for the SMA solution, 38 for the LMA solution and 15 for 
the VAC solution. 



4.2 Future 

Having performed our global analysis of the recent solar neutrino data within the frame- 
work of our MHD-RSF solution to the solar neutrino problem, we are in a position to 
calculate also the expected values of a number of observables to be measured by future 
solar neutrino experiments, such as SNO or Borexino. This task has been developed for 
the case of oscillation-type solutions to the solar neutrino problem in ref. ||33|| . Here we 



will consider our alternative MHD-RSF solution described in sections 2 and 3, because 
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Solution 


Am 2 


B_\_max 


Xlin ( Prob %) 




MED - RSF a 
MED - RSF S 


1.1 x 10~ 8 
1.1 x 1(T 8 


80 
75 


25.7 (32) 
30.5 (14) 


this work 
this work 




Am 2 


sin 2 (2$) 


xLn (Prob %) 


Ref. 


SMA a 
LMA 
LOW 

SMA S 
VAC 


5.2 x 10" 6 
2.4 x 10~ 5 
1.0 x 10~ 7 
5.2 x 10~ 6 
4.4 x 10~ 10 


4.7 x 10~ 3 

0.78 

0.93 
4.7 x 10~ 3 
0.9 


29.7 (16) 
27.0 (26) 
32.0 (10) 
32.0 (10) 
34.3 (6) 


I 
1 
1 
1 


27 
27 
27 


iz 


^ 

31 
31 

32 


1 


no-osc 






87.9 (6 x 10" 7 ) 


1 


27 


1 



Table 5: Best fit points and the corresponding probabilities for different solutions to the 
solar neutrino problem. The top row corresponds to the MHD-RSF solution presented 
here. 



of its theoretical elegance and the good quality of the global fits it provides. Again, the 



results of refs. [EyL 31, 32 1 will allow us to compare quantitatively our simplest MHD- 



RSF predictions with those associated with the vacuum (VAC) and MSW-type (SMA, 
LMA and LOW) solutions to the solar neutrino problem using the same well-calibrated 
theoretical calculation and fit procedure. 

We determine the expected solar neutrino rates at SNO using the cross sections for 



the CC and NC ud reactions given by ref. and the preliminary SNO collaboration 
estimates for the energy resolution, absolute energy scale and detection efficiencies. For 
definiteness we adopt the most optimistic threshold energy of 5 MeV which should be 



reached by the collaboration ||35|| . We perform these calculations at the best-fit points 
which we have determined in the present paper, using 90 and 99 % CL error bars. For 
definiteness we have considered the global best fit points and local minima for B± max < 
lOOkG given in table [| for the case ku = 6 and r = 0.6 and active-active MHD-RSF 
conversions. 

We have calculated the neutral-to-charged-current event ratio (NC/CC for short) and 
our results are presented in Fig. (|7]). Our predictions for the oscillation solutions agree 



relatively well with those of p3 |. The agreement is not perfect because we use the full 
zenith angle dependence in the analysis of the solar neutrino data instead of simply the 
day-night asymmetry employed in ref. f3~3| . The size of the error bars displayed in Fig. 



17 



o 
o 



7 r 



6 

5: 



LMA 



LOW 



MHD— RSF 4 



VAC 



SMA A 



SMAs 



MHD-RSF. 



Scenario 



Figure 7: Neutral-to-charged-current event ratio expected at SNO for different solutions 
to the solar neutrino problem at 90% CL and 99% CL. The no-oscillation or SM case is 
denoted by the horizontal line at one. 

(|7|) arises from the variation in the values of neutrino oscillation parameters, rather than 



from statistical and theoretical uncertainties, which are negligible |33| . 

Clearly from Fig. (]?]) we see that there is a substantial overlap between our MHD- 
RSF predictions and those found for each of the oscillation solutions (SMA, LMA, LOW, 
VAC). The overlap is especially large between the LMA and the MHD-RSF solutions. 
Taking into account the present theoretical uncertainties and a reasonable estimate of the 
experimental errors attainable, it follows that an unambiguous discrimination between 
our MHD-RSF solution and the neutrino oscillation-type solutions to the solar neutrino 
problem on the basis of the averaged event rates seems rather difficult. The expected 
features of the MHD-RSF recoil electron spectrum will be discussed elsewhere [[36 . 



5 Discussion & Conclusions 

We have re-analysed the status of resonant spin-flavour solutions to the solar neutrino 
problem in the framework of analytic solutions to the solar magneto-hydrodynamics 
equations, using the most recent global set of solar neutrino data. We have shown 
that our procedure is quite robust in the sense that the arbitrariness associated to the 
magnetic field profile has been almost eliminated due to both mathematical consistency 
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and physical requirements. Effectively our analysis substitutes neutrino mixing by a 
single parameter B± max characterizing the maximum magnitude of the magnetic field 
inside the convective region. The value of ku characterizing the maximum number of 
individual modes combined in a realistic profile and the parameter ro characterizing 
the location of the convective region are severely restricted. The allowed kyi values are 
restricted by ohmic dissipation arguments to be lower than 10 or so, and we have found 
that our solar neutrino fits are pretty stable as long as ku exceeds 5. Moreover our fits 
are pretty stable within the relevant narrow range for ro- This way we obtain effective 
two-parameter global fits of solar neutrino data for static MHD solutions characterized 
by Am 2 and B± max , since the magnetic field profile is essentially unique. This enables us 
to compare their quality with that of the fits obtained for the favored neutrino oscillation 
solutions to the solar neutrino problem. Adopting the Standard Solar Model we have 
found the MHD-RSF fits to be slightly better than the oscillation fits, though not in a 
statistically significant way. We have also analysed the prospects to distinguish our best 
MHD-RSF solution from the oscillation solutions (SMA, LMA, LOW, VAC) at future 
solar neutrino experiments. Both in the comparison of the present status of different 
solutions of the solar neutrino problem, as well as in their future predictions at SNO we 
have used a common well-calibrated theoretical procedure and statistical criteria. Taking 
into account the present theoretical uncertainties and the expected experimental errors 
attainable, an unambiguous discrimination between our MHD-RSF solution and the 
neutrino oscillation-type solutions to the solar neutrino problem at the SNO experiment 
seems rather difficult. On the other hand better measurements of rate-independent 
solar neutrino observables such as the day-night asymmetry and seasonality would be 
potentially useful, since our MHD-RSF predictions differ from the expectations of the 
oscillation schemes. For example, seasonality is expected to be smaller |36|] in our MHD- 
RSF solution than in MSW []37| or just-so oscillations j32]. On the other hand the 
day-night asymmetry of the MHD-RSF solution is negligible, in contrast with the MSW 



solutions [36 



Note, however, that the complete MHD-RSF solution is characterized also by a non- 
zero neutrino flavour mixing. This gives it the potential to be discriminated from the 
oscillation-type solutions |36| . The most distinctive signal expected in this case consists 
of solar anti-neutrinos, which would provide a clear signal in water Cerenkov exper- 



iments |(38| . Moreover, for large enough neutrino mixing one expects also a sizeable 
suppression of the rates for pp neutrinos, potentially testable at the GNO experiment. 
Last but not least, the possible time dependence of the charged current signal due to 
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solar cycles still remains as a possible tool to discriminate the MHD-model from the 
oscillation schemes. 

Note added: As we finished our paper there appeared the paper E. K. Akhmedov 
and J. Pulido, [hep-ph / 00051 73| , which also considers predictions for SNO observables 
in the RSF scheme employing the phenomenological magnetic field profiles they used 
previously in ref [|15| . 
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Angular Range 


Data; ± (Tj 


Day < cos 6 < 1 
Nl -0.2 < cos 9 < 
N2 -0.4 < cos 9 < -0.2 
N3 -0.6 < cos9 < -0.4 
N4 -0.8 < cos6 < -0.6 
N5 -1 < cos 9 < -0.8 


0.463 ±0.0115 
0.512 ±0.026 
0.471 ±0.025 
0.506 ±0.021 
0.484 ±0.023 
0.478 ±0.023 



Table 6: Super-Kamiokande Collaboration zenith angle distribution of events [El. 

A Zenith and Spectrum Data Samples and Fit Pro- 
cedures 

Here we summarize here the data used and the fit procedures adopted in this paper. 

The zenith dependence data given by the Super-Kamiokande collaboration [EJ are 
shown in table |6]. 

The recoil electron spectrum data are given as In table |7| Oi^tat is the statistical 
error, ai^ exp is the error due to correlated experimental errors, (Ji >C ai is the error due to 
the calculation of the expected spectrum, and cr itUncorr is due to uncorrelated systematic 
errors. 

In our study we use the experimental results from the Super-Kamiokande Collabo- 
ration on the recoil electron spectrum on the 18 energy bins including the results from 
the LE analysis for the 16 bins above 6.5 MeV and the results from the SLE analysis for 
the two low energy bins below 6.5 MeV, shown in table [7|. 

Notice that in table |7] we have symmetrized the errors to be included in our x 2 
analysis. We have explicitly checked that the exclusion region is very insensitive to this 
symmetrization. We define x 2 f° r the spectrum as 

-nth nth 

yl = V (a — R exp )a~ 2 (a °- R exp ) (19) 

AS 1 v" s P pBP98 l i l u n v" s P pBP98 l J > V ic > 

i,i=l,18 3 

where 

&ij = "ijK&i^tat ± ®i,uncorr) ± ~i ,expO 'j ,exp ± &i,calO ~j,cal \^) 

Again, we introduce a normalization factor a sp in order to avoid double-counting with 
the data on the total event rate which is already included in x\- Notice that in our 
definition of x% we introduce the correlations amongst the different systematic errors in 



21 



Energy bin 


flnfn -1- 

JJatai ± a^ s tat 


&i,exp l/oj 


&i,cal I/O J 


rr 

Vi,uncorr \>v) 


o.o Mev < hj e < o iviev 


n A 70 X n (1Q7 

U.4 IZ± U.Uo ( 


1 1 
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n q 
U.o 
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4.U 
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n q 
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n /i 07 J- n noo 
U.4Z / ± U.UZZ 


1 Q 
l.O 


n q 
U.o 


1 7 
1. 1 


1 MeV < Hi e < i .0 MeV 


n /i fin _i- n noo 

u.4oy ± u.uzz 


1 Q 
l.O 


n £ 
U.O 


1 7 
1. / 


/.0 MeV < -C/ e < o MeV 


n £ 1 j- n noo 
U.01O ± U.UZZ 


1.0 


n 7 
U. ( 


1 7 
1. / 


o IVleV < -C/ e < o.O MeV 


n /i sfi -i_ n noc; 
U.4oo ± U.UzO 


l.O 


n n 

u.y 


1 7 
1. f 


o.o Mev < -c/ e < y Mev 


n aaa -i_ n noc; 
U.444 ± U.UZO 


O 

z.z 


1.1 


1 7 
1. ( 


y Mev < -c/ e < y.o Mev 


n /i -i_ n noc; 
U.404 ± U.UZO 


O K 

z.o 


1.4 


1 7 
1. 1 


y . o ivie v <. .c/ e <. i u ivie v 


n e;i p. 4- n non 
u.oio it u.uzy 


O Q 

z.y 


1 7 
1. ( 


1 7 
1. ( 


10 MeV <E e < 10.5 MeV 


0.437 ±0.030 


3.3 


2.0 


1.7 


10.5 MeV < E e < 11 MeV 


0.439 ±0.032 


3.8 


2.3 


1.7 


11 MeV < £ e < 11.5 MeV 


0.476 ±0.035 


4.3 


2.6 


1.7 


11.5 MeV < E e < 12 MeV 


0.481 ±0.039 


4.8 


3.0 


1.7 


12. MeV < E e < 12.5 MeV 


0.499 ± 0.044 


5.3 


3.4 


1.7 


12.5 MeV < E e < 13 MeV 


0.538 ±0.054 


6.0 


3.8 


1.7 


13 MeV <E e < 13.5 MeV 


0.530 ±0.069 


6.6 


4.3 


1.7 


13.5 MeV < E e < 14 MeV 


0.689 ±0.092 


7.3 


4.7 


1.7 


14 MeV < £ e < 20 MeV 


0.612 ±0.077 


9.2 


5.8 


1.7 



Table 7: Recoil energy spectrum of solar neutrinos from the 825-day Super-Kamiokande 
Collaboration data sample H]. 
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the form of a non-diagonal error matrix in analogy to our previous analysis of the total 
rates. These correlations take into account the systematic uncertainties related to the 
absolute energy scale and energy resolution. 

The general expression of the expected rate in the presence of oscillations R th in a 
bin, is given from Eq.fllBD but integrating within the corresponding electron recoil energy 
bin and taking into account that the finite energy resolution implies that the measured 
kinetic energy T of the scattered electron is distributed around the true kinetic energy 
T' according to a resolution function Res(T, T') of the form 



Res(T, T) 



2ns 



cxp 



(T-T 
2? 



A2 



where 



s = s JT'/MeV 



(21) 



:22) 



and So = 0.47 MeV for Super-Kamiokande j|, |39|. On the other hand, the distribution 
of the true kinetic energy T' for an interacting neutrino of energy E v is dictated by the 



differential cross section da a (E„, T')/dT', that we take from [P9| . The kinematic limits 
are 



< T' < T\E U ) , T\E V ) 



E u 



(23) 



1 + m e /2E u ' 

For assigned values of s , T min , and T max , the corrected cross section a a (E) (a = e, x) is 
given as 



dT 



T(& V fl e 3 (T,T0^- r) 



dT 



(24) 
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